The Computational Power of Symmetric Hamiltonians 
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The presence of symmetries, be they discrete or continuous, in a physical system typically leads 
to a reduction in the problem to be solved. Here we report that neither translational invariance nor 
rotational invariance reduce the computational complexity of simulating Hamiltonian dynamics; the 
problem is still BQP complete, and is believed to be hard on a classical computer. This is achieved 
by designing a system to implement a Universal Quantum Interface, a device which enables control 
of an entire computation through the control of a fixed number of spins, and using it as a building- 
block to entirely remove the need for control, except in the system initialisation. Finally, it is shown 
that cooling such Hamiltonians to their ground states in the presence of random magnetic fields 
solves a QMA-complete problem. 



I. INTRODUCTION 

Within the study of quantum computation, it is desir- 
able to find 'natural' problems for a quantum computer 
to solve i.e. those that exhibit intrinsically quantum prop- 
erties which elucidate the power of the device. For ex- 
ample, Feynman first proposed the quantum computer as 
a device that naturally simulates Hamiltonian dynamics. 
More recently, it has become apparent that finding the 
ground state energies of Hamiltonians is QMA-complete 
IE S @, EL HI i-e. this is a natural problem for the class 
where solutions can be efficiently verified on a quantum 
computer. The question that naturally arises with re- 
gards to both of these problems is that of the minimal 
properties that any such Hamiltonian must posses. The 
same question, viewed from another perspective demands 
when we should expect efficient classical approximations 
to the properties and dynamics of Hamiltonians. 

It is often found that the presence of symmetries, such 
as translational or rotational invariance, vastly simplify a 
problem. It may be hoped that accounting for these prop- 
erties could make classical simulations of ground state, 
thermal or dynamic properties tractable and yet perti- 
nent to physical phenomena. A straightforward exam- 
ple of this reduction is the description of a two-qubit 
mixed state p. To do this ordinarily requires 15 real pa- 
rameters. However, if the state is rotationally invariant, 
i.e. {U ®U)p(U^ ®W) = p for all single-qubit unitaries U, 
then the family of possible states, the Werner states @, 
is parametrised by a single number. Similarly, under the 
discrete permutation symmetry, SWAP • p ■ SWAP = p, 
the family of states only contains 10 real parameters. The 
hope that the introduction of translational invariance re- 
duces the complexity of finding ground states underpins 
studies such as Matrix Product States @, B Q • 

In this paper, we are primarily concerned with the ef- 
fects that these symmetries have on the ability to effi- 
ciently simulate Hamiltonian dynamics, and report that, 
in fact, the symmetries have no bearing; the problem is 
BQP-complete i.e. as hard as any quantum computation 
is to simulate on a classical computer. We will show this 
by constructing Hamiltonians that implement arbitrary 



quantum computations. 

The first steps towards incorporating translational in- 
variance, for both Hamiltonian evolution and ground 
state properties, were taken in Q, which traded the 
spatial variation for another property such as the local 
Hilbert space dimension, which grew as poly(A) for N 
spins. Necessary to this construction was the inclusion 
of both a time label (clock) and position label at each 
site, so that all the information was available locally. 
Processing was then achieved by implementing a read 
head that moved backwards and forwards over the sys- 
tem data. Here, we use global control (GC) schemes 
[ToL [Il| to remove the necessity of both time and po- 
sition labels. In Sec. [Til we w iU describe how to apply 
global commands through Hamiltonian evolution leading, 
in Sec. IIIIl to the realisation of a translational invariant, 
nearest-neighbor Hamiltonian of fixed local Hilbert space 
dimension that implements arbitrary quantum computa- 
tions, thus implying the classical intractability of simu- 
lation of the dynamics, of which there is already some 
evidence [13]. The construction is easy to motivate be- 
cause a GC scheme typically works by repeated applica- 
tion of a finite set of pulses "A" , "B" etc., which are local 
gates applied uniformly to all qubits in the system. As 
such, we immediately lose the need for spatial resolution 
in our Hamiltonian. To remove the time resolution, the 
program sequence is written on the initial state of some 
of the spins instead of encoding it in the Hamiltonian. 
Proceeding to a proof of QMA-completeness for ground 
state properties still requires the introduction of spatially 
varying magnetic fields (Sec. IIV|). but this in turn has se- 
vere implications for the cooling of physical systems in 
the presence of random external fields. 



II. A UNIVERSAL QUANTUM INTERFACE 

Lloyd and co-workers proposed the concept of a Uni- 
versal Quantum Interface (UQI) [la ]. They proved that 
through manipulation of a single spin which is coupled 
to a larger system, the entire quantum system can be 
controlled, such that a quantum computation can be im- 



plemented on it. The proof was non-constructive, and 
there has recently been some interest in how such a de- 
vice might be implemented An example of one of 
the potential benefits of such a scheme is that it would 
in principle allow one to isolate the bulk system from the 
environment so that it's much less susceptible to noise. 
In comparison to schemes for Hamiltonian evolution, the 
UQI protocol also incorporates the ability to prepare the 
initial state of the system, although we will not explore 
that aspect here. Control theory proves the existence of 
solutions and offers numerical techniques to determine 
high-fidelity control 14]. In contrast, we realise an exact 
UQI constructively. This will allow us to introduce most 
of the concepts required throughout the paper. Indeed, 
we will use the UQI as a building block in the design of 
our computing Hamiltonians. 

Consider a linear chain of N 4-dimensional spins. 
These 4-level spins can be assigned a structure as the 
tensor product of two qubits, labelled a and q. Systems 
q will contain computational qubits, and systems a will 
be used as a 'read head'. If the spins are coupled by an 
interaction 
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where Sl i+1 is the swap operation between computa- 
tional qubits i and i + 1, and X and Y are the stan- 
dard Pauli operators, then control of just spins 1, 2 and 
N is sufficient to realise universal quantum computation. 
We can understand this by transforming the Hamiltonian 
into a state transfer system. We are interested in the sub- 
space when all the a qubits are in the state |0), except for 
qubit i which is in the state 1 1) . We choose to denote that 
as \i) a \4>i) q where \4>i) q is the state of the computational 
qubits in that step, and satisfies \ipi) — S^-i^ 
With this definition, one can see that 

ffuqi |i>° \^) q = \i - l) a \A-x) q + \i + l) a \^+i) q 

(|0) a and \N + l) a are taken to be 0) which is just the 
same as the model of state transfer studied by Bose [l5[ , 
except that one needs to apply local magnetic fields on 
spins 1 and N. In particular, Bose studied the evolution 
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W) a \"4>n} 9 [IE EH, where one can achieve 



an arrival probability of 0(N~ 2 / 3 ) within a time 0(N) 
[To] ]. Moreover, this transfer is heralded - by measuring 
system a on qubit N, we know whether or not the trans- 
fer has occurred without disturbing the computational 
qubits. If the transfer is not finished, we simply wait and 
try again. With this strategy, one can achieve an evo- 
lution which fails with probability less than s in a time 
0(M 5 / 3 log(l/e)) []]|. The setup is depicted in Fig.[TJ 

The computational scheme now proceeds as follows. 
We assume the state to be initialised with all the a qubits 
in |0), and the computational qubits are in state \ipi). 
The read head is initialised by placing the a qubit of 
spin 1 in the |1) state, and we perform a heralded trans- 
fer to spin N, which means that when it is observed 
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FIG. 1: The basic scenario of a universal quantum interface. 
At each site, there are 2 qubit systems, a and q. All the as are 
initialised as |0). A |l) a is input at site 1, and its arrival at 
site N is monitored. The hopping Hamiltonian is such that as 
the excitation moves between site i and i + 1, a unitary L/^j+i 
is enacted on the computational qubits q. Thus, upon arrival 
at site N, the unitary Un-i,n ■ ■ ■ U\,2 has been implemented. 
Furthermore, this transfer is efficient. 



that the excitation has arrived at \N) a , we reset the spin 
to |0), and the computational qubits have changed to 
Sn-\,n ■ ■ ■ <S , 2,3<5'i,2 which is simply a cyclic permu- 

tation of the computational qubits. This can be repeated, 
allowing any computational qubit to be placed on spin 1 . 
Similarly, by starting the excitation on spin 2, and re- 
moving it from spin N, the computational qubits 2 to N 
undergo a cyclic permutation, allowing us to place any 
single qubit on spin 2 without disturbing spin 1. One we 
have any arbitrary pair of qubits on spins 1 and 2, which 
requires no more than cycles, our control of these two 
qubits allows the implementation of an arbitrary one- or 
two-qubit gate, which is sufficient for universal quantum 
computation. 



A. A Toolbox of UQI Gadgets 

Having demonstrated a relatively simple way to imple- 
ment a UQI, we will now show how to modify the scheme 
to improve its capabilities. Our first observation is that 
we can replace the swap gate, Sf i+1 in Eqn. (fTJ) with 
any unitary C/ 2 9 i+1 , and the transfer |1)° — * |1)^ imple- 
ments the operation Un-i,n ■ ■ ■ ^2,3^1,2 The only 
difference is that we must ensure that the Hamiltonian 
remains Hcrmitian, which we do by writing it as 

ffffii = E( ff rV+i) ® u q l+1 + K + v-) ® uj q +v 



The Hamiltonian remains translationally invariant pro- 
vided Ui t i+i is the same as U\^, just acting on different 
qubits. 

The next step is to see how to implement L^^+i only 
on every second qubit. We achieve this by increasing the 
dimension of system a to 3. The idea is that instead of 
having a Hamiltonian where the read-head hops to the 
right, implementing U as it goes, it will alternate its value 
between 1 and 2 as it hops, and will only implement U if 
it's hopping from 1 to 2, and not 2 to 1. Thus, the new 
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Hamiltonian reads 

4qi = E l° 2 ) (!0|m + i®^ + i + |01) (20|? i+1 ®l < * i+1 +fc.c. 

i 

When we initialise the read- head on spin 1, whether it's 
in the state |1) or \2) determines which pairs of compu- 
tational qubits the U is applied to, pairs (2i — l,2i) or 
(2i, 2i+l), and also what state we should be checking for 
on spin N. If N is odd, the read-head undergoes an even 
number of steps and arrives in the same state it started 
on spin 1. 

Finally, we might like the choice of applying more than 
one different U, say U and V. Again, we can achieve this 
with the system a having 3 levels. 

H v% = E I 01 ) aCm^+i+lOS) (20\i,i + i®V? i+1 +h.c. 

i 

By initialising the read-head on spin 1 as |l) a , the |l) a 
hops from site to site, and as it hops it implements U. 
Similarly, initialising spin 1 in the state |2) causes V to 
be implemented. 

These constructions are readily combined so that if 
we want to implement, say, two different global com- 
mands, each of which only acts on every second qubit, 
then we can construct a translationally invariant Hamil- 
tonian that does it, using a read-head of dimension 5, and 
thus an overall spin dimension of 10. It is notationally 
convenient at this stage to decompose the 5-dimcnsional 
read-head at each site i into 

(aj ® n) © Hi. 

The state |n) (the equivalent of |0) a in the previous no- 
tation), is the state which all the read- head systems are 
initialised in, except for one, and is used to indicate the 
absence of the read head. The systems a« and are both 
qubits. <Zj indicates if the read-head is 'active', and rj 
contains the 'program information', i.e. which of the two 
gates to implement. When we propagate the read-head 
to the right, the interaction that is implemented on the 
data qubits is conditioned on the value of the read-head, 
and whether the read-head is active. As the read-head 
hops, we flip the active setting, so that it only applies an 
operation on every second qubit. 

III. COMPUTATION BY HAMILTONIAN 
EVOLUTION 

In [l7|, a global control scheme was developed based 
on just two nearest-neighbor gates, SWAP (S^j+i) and 

G M+ i = 1 [Z - iY)/V2. 

These gates need to be applied to distinct qubit pairings 
(2i — 1, 2i) and (2i, 2i + 1) across the entire lattice. Evi- 
dently, this ties in very usefully with our UQI construc- 
tion. If we implement these two gates within the Hamil- 
tonian, then the 10 level system is capable of universal 
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FIG. 2: Schematic of the Hamiltonian's mechanism. At each 
site, there is a label to specify if that site is a computational 
qubit or program qutrit. One location in the program bus is 
marked as the active region, and that value is stored in the 
read- head qubit r%. The read-head applies a unitary to each 
computational qubit consecutively, controlled off its value. 
When it reaches the active program region, it exchanges its 
data with the next step in the program. 

quantum computation, where we only need single-spin 
control of spins 1 and N to perform the transfer of the 
read-head, and to specify the program sequence (i.e. the 
order in which global commands are implemented). 

Our aim is now to construct a Hamiltonian that im- 
plements the entire evolution of the computation without 
even this basic level of control, but instead retaining the 
ability to prepare the system in some initial state. The 
first step is to remove the asymmetry between operations 
at either end of the chain; we introduce periodic bound- 
ary conditions and a special marker state. The read-head 
will be emitted from one side of the marker state, and 
will arrive at the other side. It will then be down to 
our preparation of the initial state to select where this 
marker state is. The second step is to incorporate the 
program sequence. Again, we will do this by writing 
the additional information in the initial state. This will 
require an additional two-level system at each site to de- 
note whether that site constitutes a computational qubit 
or a spin that holds the program data. The previously 
mentioned marker will then be used to denote which of 
the program spins is the one that's being actively imple- 
mented. Again, the read-head will be emitted from one 
side of this marker, but when it arrives on the other side, 
it will move the marker onto the next program spin (see 
Fig. The natural start and end points to the compu- 
tation are when the active program marker is at cither 
end of the program bus, and are detected by a change in 
the li label between program and data spins. 

To be specific, consider a ID chain of spins of local di- 
mension 31, which can be decomposed into several sub- 
systems 

Si = (qi <S> k ® ((oi ® n) © m)) © TOj. (2) 

The system qi is a qutrit, serving two different purposes 
depending on the label of the qubit system k . If the label 
is |0), then |0) 9 , |1) 9 encode a computational qubit (the 
system q in the UQI construction), otherwise the qutrit 
qt contains program information - "skip" , G or S. There 
is a single state which is not used yet, \2) q |0)', and is 
reserved as the marker to denote which program trit is 
active. The single level is used to help moving the 
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marker over a 'skip' label. 

Whether a particular global gate works on pairings of 
qubits (2i — 1, 2i) or (2i, 2% + 1) is solely determined by 
the alignment of the active program trit with respect to 
the start of the block of data qubits, which is why we 
require skip; such that the relative alignment changes. 
Assuming that there are an odd number of spins in the 
system, if a read-head leaves a program trit in the active 
state, it returns in the inactive state. Consequently, an 
inactive read-head in the location of the active program 
trit can be used to indicate that the read-head should 
move the active trit to the next one. 

The main term in the Hamiltonian H prop involves read- 
head propagation, and comes directly from the UQI con- 
struction. As such, we will not repeat it here. We must 
additionally incorporate a term to stop the read-head 
propagating if the state \2) q |0)' is present i.e. 

H; rop = H prop (t - |20) (20|f )+ 
J2 X \n20) (zl20|[ O9/ ® (\x0) ra (n| ® l ql ) l+ i + h.c. 

Next, program manipulation (when the read- head is in 
the neighborhood of the currently active program state): 

H; rog = ^(|n)(x0|™®|x + l)(2|^|l)(0|') l 

x,y 

®(\yl) ra (n| ® |2) (y + l\ q ® |0) + h.c. 

except that this doesn't (yet) handle the skip. First, if 
we're in a region where we've just arrived back from doing 
a loop, and should be moving onto a skip label 

Hli = ^(|n)(x0|™®|x + l)(2| 9 ®|l)(0| / ) l 

X 

®(|m) (n\(0l\ ql ) l+1 + h.c. 
and, secondly, the step over the skip label 

H l S 2 = £(l n >l 01 > a, M)«® 

X 

(\xl) ra (n\ ® \2) (x + l\ q ® |0) + h.c. 

Finally, the total Hamiltonian is 

Ht = H prop + H prog + H l sl + H l s2l (3) 

i 

which is entirely a sum of two-body terms, so it can be 
represented as J2i Moreover, for the permutation 

operator 

30 

P= ^ 1*1*2 ■ ■ -In) (*2«3 ■ ■ • *JV«lj , 


we have that 

PH T P ] = H T 



i.e. the Hamiltonian is translationally invariant. 

If there are N computational qubits in the system, then 
any efficient quantum algorithm is described by poly(TV) 
bits, and the total number of qubits in the system is 
M, the sum of these two. Again, we can invoke the 
fact [HI that after an evolution time O(M), the prob- 
ability of having successfully completed the computation 
of 0(M" 2 / 3 ), which is an efficient implementation. So, 
this translationally invariant, nearest neighbor Hamilto- 
nian evolution can implement any quantum computation, 
starting from a separable state. The fact that it starts 
from a separable (although not translationally invariant) 
state is important since it ensures that we are not encap- 
sulating the difficulty of the problem within the prepa- 
ration of the initial state. We conclude that the problem 
is BQP-hard. However, it is also known how to simulate 
Hamiltonian evolution on a quantum computer [l8l |. so 
the problem is BQP-complete. Thus, as strongly as we 
believe that quantum computation is more powerful than 
classical computation is how strongly we believe that sim- 
ulation of Hamiltonian dynamics, even under the trans- 
lational invariant restriction, is hard to simulate on a 
classical computer. 

A. Rotational Invariance 

So far, we have seen how a ID Hamiltonian with fixed 
local Hilbert space dimension and fixed range interac- 
tions, which is translationally invariant, can implement 
an arbitrary quantum computation. However, this dis- 
crete symmetry is not nearly as restrictive as the contin- 
uous symmetry of rotational invariance, which requires 
the Hamiltonian to satisfy 

for all single qubit unitaries U. We will now show how to 
build this into the Hamiltonian, retaining translational 
invariance and the ability to perform arbitrary quantum 
computations. The first step is to take the 31 dimen- 
sional construction Ht, Eqn. ([3]), and replace each spin 
with 10 qubits. Between these 10 qubits, there are several 
decoherence-free subsystems [l9[ . For an N qubit system 
(iV even), there are decoherence-free subsystems which 
have (^ 2 ) (2j + l)/(N/2 levels for j = 0... N/2. 

These subsystems enable the storage of quantum infor- 
mation in a way that is not affected by collective deco- 
herence U® N . Thus, encoding within one of these sub- 
systems stores the information in a rotationally invariant 
way. We select any one of the 4 subsystems (N = 10) 
that is large enough to encode the 31 levels. Transcribing 
Ht into this new form automatically makes it rotation- 
ally invariant, although instead of being translationally 
invariant, it is periodic, with a repetition length of 10 
qubits. We can thus write it as 

M-l 

Hr= fyl0(i-l)...10i-l),(10i...l0(i+l)-l); 
i=l 
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(c) (f i s f r r~sn F r~sn f T~sn ^ r~s~r n r~5~) 

FIG. 3: (a) Schematic depiction of a ID array of 31 di- 
mensional spins (s). The black regions denote the nearest- 
neighbor interactions of the Hamiltonian. (b) The same sys- 
tem made rotationally invariant by encoding the states of the 
spin in logical states on a block of qubits. The logical spins, 
and the periodic Hamiltonian interactions are depicted. To 
make the system translationally invariant requires the inclu- 
sion of terms (gray) where the interactions do not align with 
the logical qubits of the state, (c) To regain translational 
invariance in the Hamiltonian, we introduce a flag state (f) 
before each block of logical qubits, denoting the start of that 
block. 



denoting the blocks of logical spins. A summary of this 
is depicted schematically in Fig. [3] Note that the ini- 
tial state that the Hamiltonian evolution acts on must 
also be encoded. However, it is encoded into fixed sized 
blocks which are separable from each other. Thus, the 
initial state can still be efficiently represented on a clas- 
sical computer, so we have not transformed the problem 
of simulation into the preparation of the initial state - 
it's still contained within the Hamiltonian evolution. 

In order to reintroduce translational invariance, we 
want to incorporate a local patterning of states that en- 
ables us to detect the alignment of the blocks of qubits 
(Fig. 02(c)). The technique that we use is much clearer 
if we concentrate, initially, on restricting the arbitrary 
rotations to rotations about a single axis, U = e . Our 
problem is that a translationally invariant Hamiltonian 
will be made up of sums of terms, each of which acts 
on a block of qubits, comprised of two logical spins. We 
need to make sure that if a Hamiltonian term is not per- 
fectly aligned with the block- wise patterning of the initial 
state, then it does not contribute to the evolution. To do 
this, we introduce a patterning of the qubits which is 
still rotationally invariant, and yet flags the start of each 
block of spins. For Z-rotation invariant states, this can 
be done by taking each set of qubits that constitutes a 
logical spin, and introducing a qubit in the |0) state be- 
tween each of them. This means that in the initial state, 
there is never a pair of neighboring qubits in the |11) 
state. Thus, we use this to flag that a block of logical 
qubits is starting. So, each logical spin now constitutes 
23 qubits, of the form \HOq1Oq2Oq3Oq4O . . . qioO) where 
\liQ2 ■ ■ ■ Qw) wa s the previous logical spin encoded into 
a decoherence-free subsystem, and the Hamiltonian is of 
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FIG. 4: The rotationally invariant state that flags, and the 
Hamiltonian that detects the flag, at the start of the block 
of 10 qubits encoding a logical spin within a decoherence free 
subsystem. Section A makes sure that when the Hamiltonian 
and herald state are offset by a single qubit, the overlap is 0. 
Section B, where the states P3 are repeated 8 times, handles 
a relative shift of an even number of qubits. Section C does 
the same for an odd number of qubits, and the states P3 are 
repeated 9 times. The difference in the number of repetitions 
of P3 handles an edge effect that arises otherwise. 



the form |11) (11| 1)2 /i(4,6,8...22), (27,29.. .45) ■ 

For Z-rotations, we have the advantage that the states 
|0) and |1) are rotationally invariant (but not superpo- 
sitions of them). For arbitrary rotations U, the con- 
struction of suitable flag states is much more involved, 
and is based on the fact that \ip~) = (|01) - |1Q))/a/2 
is a rotationally invariant two-qubit state. Hamilto- 
nian terms can be constructed which detect the pres- 
ence, P 1 = or absence, P 3 = 1 — \if)~) {tl>~\, 
of such a state. There will be a correspondence between 
the Hamiltonian term that detects the flag state, and the 
flag state itself, the only difference is that if the Hamilto- 
nian includes a term P 3 , the flag state can be made out 
of any two-qubit pure state orthogonal to such as 
1 00). We proceed by observing that it is relatively easy 
to suppress misalignments between the flag state and the 
Hamiltonian term when the misalignment is by an even 
number of qubits, one simply ensures that a in the 
flag state and a P 3 in the Hamiltonian align. This trick 
can be repeated for a misalignment by an odd number of 
qubits greater than 1. There are two concerns remain- 
ing. Firstly, whether there any edge effects arising and, 
secondly, how to deal with an offset of just one qubit. 
The first concern is overcome simply by using a large 
enough flag state, which is larger than the 10 qubits in 
the logical spin. An offset of one qubit is handled by in- 
corporating four additional qubits in the flag state in the 
form -P/ 3 <8 -R 3 4- The interleaving of the two projectors 
ensures overlap after just a single shift. The entire state 
is depicted in Fig. 2] and requires 43 qubits in total. The 
projector onto the state is written as 

7 

Pl fla .43 = A>A>^>0^ + 1,W® 

i=0 

8 

®P 2 3 ® 134 ® <S> P 2i+25,2i+26- 
i=0 

Thus, the overall translationally and rotationally invari- 
ant Hamiltonian acts on blocks of 106 qubits which are 
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local on a ID lattice. It is of the form 
ftf-i 



H 



Epfiag 
^53(i-l). ..53(i-l) 



P r 



flag 



RT ~ / , 1 S3(t-1)...S3(i-1)+42 ^ J 53i...53i+42 

i=l 

®^(53(i-l)+43...53i-l),(53i+43...53(i+l)-l)) 

where the hnn are the same as in Hr. The majority of 
the cost in terms of the range of the Hamiltonian terms is 
due to the flag state, which we have made little attempt 
to optimise; the important element is that the range of 
the terms is independent of N. 



IV. GROUND STATE PROPERTIES IN THE 
PRESENCE OF MAGNETIC FIELDS 

In previous studies, Hamiltonian evolution has been 
used as a basis for classifying the problem of finding 
ground state energies of Hamiltonians with similar prop- 
erties as QMA-complete. In the present case, this is not 
expected to be possible as there seems to be no way to 
encode the verifier's computation while retaining transla- 
tional invariance. However, by breaking the translational 
invariance of Ht, one arrives at similar results to j4|, but 
only using local magnetic fields. To achieve this, we need 
to add several energy penalties; to detect the solution to 
the verifier circuit, to initialise ancillas in |0) and to pre- 
pare the initial state of the program tape for a specific 
computation corresponding to the verifier of the QMA 
problem. To implement these penalties, we need to be 
able to locally detect that we are either at the beginning 
or end of a computation, requiring an increase in the lo- 
cal Hilbert space dimension, such that the system n has 
2 levels. |0) n will be used as before, to indicate that 
the read-head is inactive. |1) is a program command 
that will only get used once, as the first command. It is 
not necessary to program it, because two-body terms can 
readily detect the transition between the data and pro- 
gram spins. A further change is that the system a must 
be increased to dimension 3, leaving the overall Hilbert 
space dimension as 49. The extra level in a serves a dual 
purpose. Firstly, it can be used in the same way as |l) n , 
but to indicate the end of the computation, such that we 
can penalise the output qubit. Secondly, it is used to help 
ensure that the correct computational sequence occurs. 
We will adapt the Hamiltonian propagation such that if 
a read-head in either |0) a or |l) a arrives at the region 
of transition from data to program spins, it is converted 
into 1 2) , which will continue to propagate through the 
program region until it gets to the active program spin, 
where it releases its information and gets reinitialised in 
|l) a . If the read- head reaches the end of the program 
region in the \2) a state, it is deactivated, and the compu- 
tation ends. In particular, this means that if the system 
were to be initialised without an active program label, 
the computation is much shorter than it would other- 
wise have been. The Hamiltonian is readily revised to 
take these alterations into account. An energy penalty 



for when the read-head passes a particular qubit in ei- 
ther |l) n or |2) a behaves exactly like the initial and final 
penalties that we require, so we simply use penalties 



E 



i) (ii 



where the tape value should be Xi (xi = for ancillas), 
and hence Xi implies a sum over all other possible pro- 
gram (data) states, including the active label |0) \2} 9 . 
The final result term (to test the output of the verifier 
on qubit o) is similar: 



Hn 



|2) (2|>|0) (0|«. 



We have to be sure that the computation is initialised 
correctly, with all the spins correctly arranged. This is 
achieved by adding a constant term Hb- On program 
spins, this term is |0) (0|* <g) (|0) (0| + |1) (1|) 9 , ensuring 
that they are never data qubits. On data spins, this term 
is the opposite, |1) (1| ( <g> l 9 + |0) (0|'<g)|2) (2| 9 . Taking all 
of this into account, one can directly apply the projection 
lemma of Q, 

Lemma 1. Let H = Hi + H 2 be the sum of two Hamil- 
tonians operating on some Hilbert space Ti. = S + S . 
The Hamiltonian H2 is such that S is the ground state 
eigenspace (with eigenvalue 0) and the eigenvectors in S 1 - 
have eigenvalue at least J > 2||J?i||. Then, 



is 



J- 2 Lff: 



<X(H). 



X(H) denotes the smallest eigenvalue of H. 

For example, with J = 8||i?i|| 2 + 2||i?i||, one obtains 
X(Hi\s) — h A (if). We start with the total Hamiltonian 



H, 



total 



-JqHt + JbHb + JinHin + JoutHov 



Kl, 



where k = —JoX(-Ht)- Ht takes a slightly different 
form to standard proofs since instead of mapping to a 
Heisenberg chain, it maps to an XX model, and hence the 
eigenvalues are 2 cos(7rm/(M + 2)) for m — 1 . . . M + 1. 
There is an energy gap for computations that are fewer 
than M steps, as well as a gap to the first excited state, 



2 cos 
-2 cos 



M + l 

2tt 
M + 2 



2 cos 



2 cos 



7T \ C . „ 

MT-2) * I^= AE 



M + 2 



> AE 



for some constant c > 0. In the case where 'yes' so- 
lutions exist, \{H total) — JqX{—Ht) + n = 0. In the 
case where there are only 'no' solutions, we assign H 2 = 
— JqHt + JbHb + k1 to find that J > min ( Jf,, JoAE) , 
and, furthermore, \{H total ) > X(Hi\ So ) — 1/8, provided 
J > 8(J out + J in ) 2 + 2(J out + J in ), imposing a polyno- 
mial relation between Jq, Jb and J,- n , J ou t. Repeating the 



7 



process on \(Hi\ So ) with H[ = J utH 0Ut \s a , shows that 
provided J m > 8J* ut + 2J out , 

\(tt I w Jout{l-e) 2 / (M + l)7r \ 1 

where the verifier circuit of our QMA problem accepts 
the result with probability less than e. Thus, by select- 
ing J out = (M + l)sin~ 2 (j^) < c'M 3 , all the terms 
^o«t, Jim Jb and Jo are polynomial in M , and 

3 

H^total) > J ~~ £ ' 

Distinguishing the ground state energy of this Hamilto- 
nian to within l/poly(M) determines the existence of 
'yes' solutions, and thus finding the ground state en- 
ergy is QMA-complete. In comparison to Q , all of the 
spatially varying terms are local magnetic fields. An 
identical proof holds using Hut , our qubit Hamiltonian 
which is both translationally and rotationally invariant, 
although the penalties are no longer local magnetic fields 
and are, instead, 53-body. 

V. CONCLUSIONS 

Making use the GC scheme introduced in [17] . we 
have developed three main results. Firstly, the Universal 
Quantum Interface was described, and used as a build- 
ing block for the second part, which showed that the 
evolution of a fixed Hamiltonian which is translationally 
invariant on a nearest-neighbor chain and has fixed spin 
dimension can simulate any arbitrary quantum computa- 
tion, thereby suggesting that the evolution is hard to sim- 
ulate classically because it is a BQP-complete problem. 
Even simulations over short time scales, 0(AE^ 1 ^ 4 ), re- 
veal the solution since we can project onto the heralded 
outcome. We have also extended this result to include 
qubit Hamiltonians which are translationally and rota- 
tionally invariant, and still act on 0(l)-nearest neighbors 
which are local on a ID lattice. 



Finally, finding the ground state of a translationally 
invariant Hamiltonian in the presence of a specific se- 
quence of local magnetic fields is QMA-complete. This 
has a more useful interpretation as a specific example 
of a random magnetic field i.e. finding the ground state 
of a translationally invariant Hamiltonian on M spins 
in the presence of a random local magnetic field of size 
0(1/ M 2 ) is QMA-complete. This has some important 
consequences for physical scenarios involving cooling. For 
example, if one were to couple a refrigerator to a quan- 
tum system in an attempt to cool it, if the coupling is too 
strong, it would take a prohibitively long time to reach 
the ground state of the system, and that ground state 
is not the same ground state when it's not coupled to 
the refrigerator [24J. Conversely, the weaker one couples 
the system to the refrigerator, the longer it takes to cool. 
This helps to provide a motivation for the use of topo- 
logical and self-correcting systems, which are carefully 
designed such that local magnetic fields cannot affect the 
system degeneracy r .'2'.j . 

Since completing this work, we have been made aware 
of related work on evolution by translationally invari- 
ant Hamiltonians [2l[ and subsequent work (22j which 
reduces the local Hilbert space dimension in that case. 
This could presumably be used to reduce the number 
of levels required for the other results presented in this 
paper. 

In the future, we intend to examine whether the 
present work enables any useful insights into the problem 
of determining ground states for translationally invari- 
ant systems. It certainly leads to some natural conjec- 
tures which we are working to prove. Also, the technique 
for converting translationally invariant Hamiltonians into 
translationally and rotationally invariant Hamiltonians 
may be usefully applied to show that the two ground 
state problems have the same computational complex- 
ity - it is certainly true that the eigenstates of Ht map 
through to Hut- It remains to prove that the ground 
state of Ht maps to the ground state of Hjit- 

AK would like to thank F. Verstraete for useful conver- 
sions and Clare College, Cambridge for financial support. 



[1] 



[2] 
[3] 

W 
[5] 
[6] 
[7] 
[8] 
[9] 

[10] 



A. Y. Kitaev, A. H. Shen, and M. N. Vyalyi, Classical and 

quantum computation, Graduate studies in mathematics 

(AMS, Providence, Rhodes Island, 2002). 

R. Oliveira and B. M. Terhal (2005), quant-ph/0504050. 

J. Kempe, A. Kitaev, and O. Regev, SIAM Journal of 

Computing 35, 1070 (2006). 

D. Aharonov et al. FOCS07, 373 (2007). 

A. Kay, Phys. Rev. A 76, 030307(R) (2007). 

R. F. Werner, Phys. Rev. A 40, 4277 (1989). 

S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

S. R. White, Phys. Rev. B 48, 10345 (1993). 

F. Verstraete, D. Porras, and J. I. Cirac, Phys. Rev. Lett. 

93, 227205 (2004). 

S. Lloyd, Science 261, 1569 (1993). 



[11] A. Kay, PhD Thesis, http: //cam. qubit . org/users/Alastair/thesii 
(2006). 

[12] N. Schuch et al, arXiv: 0801.2078. 

[13] S. Lloyd, A. J. Landahl, and J. E. Slotine, Phys. Rev. A 

69, 012305 (2004). 
[14] D. Burgarth and V. Giovannetti (2007), arXiv: 

0710.0302; S. G. Schirmer, I. C. H. Pullen and 

P. J. Pemberton-Ross (2008), arXiv: 0801.0721. 
[15] S. Bose, Phys. Rev. Lett. 91, 207901 (2003); D. Burgarth 

and S. Bose, Phys. Rev. A 71, 052315 (2005). 
[16] M. Christandl et al, Phys. Rev. A 71, 032312 (2005). 
[17] D. Shepherd, T. Franz, and R. Werner, 

Phys. Rev. Lett. 97, 020502 (2006). 
[18] J. L. Dodd, M. A. Nielsen, M. J. Bremner and R. Thew, 



8 



Phys. Rev. A 65, 040301(R) (2002). 
[19] S. D. Bartlett, T. Rudolph and R. W. Spekkens, Rev. 

Mod. Phys 79, 555 (2007). 
[20] R. Augusiak and J. Stasinska, Phys. Lett. A 363 182 

(2007). 

[21] K. G. H. Vollbrecht and J. I. Cirac larXiv:0704.3432l 

(2007). 

[22] D. Nagaj and P. Wocian. la7Xiv:080 2.0886 (2008); B. A. 



Chase and A. J. Landahl, arXiv:0802.1207 (2008). 

[23] E. Dennis, A. Kitaev, A. Landahl and J. Preskill, 
J. Math. Phys. 43, 4452 (2002). 

[24] Although in the particular example that we have con- 
structed, the ground state of the coupled system has sig- 
nificant overlap with the degenerate ground state space 
of the original Hamiltonian. 



